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Abstract 

In [1] it has been shown that the Cabibbo angle 9c might arise from a dihedral flavor symmetry 
which is broken to different (directions of) subgroups in the up and the down quark sector. 
This leads to a prediction of 6c in terms of group theoretical quantities only, i.e. the index 
n of the dihedral group D„, the index j of the fermion representation 2j and the preserved 
subgroups indicated by ruu and ma- Here we construct a low energy model which incorporates 
this idea. The gauge group is the one of the Standard Model and D-j x z'f^^'^ serves as flavor 
symmetry. The additional jg necessary in order to maintain two sets of Higgs fields, 

one which couples only to up quarks and another one coupling only to down quarks. We 
assume that D-; is broken spontaneously at the electroweak scale by vacuum expectation 
values of SU{2)l doublet Higgs fields. The quark masses and mixing parameters can be 
accommodated well. Furthermore, the potential of the Higgs fields is studied numerically 
in order to show that the required configuration of the vacuum expectation values can be 
achieved. We also comment on more minimalist models which explain the Cabibbo angle in 
terms of group theoretical quantities, while and vanish at leading order. Finally, we 
perform a detailed numerical study of the lepton mixing matrix Vmns m which one of its 
elements is entirely determined by the group theory of a dihedral symmetry. Thereby, we 
show that nearly tri-bi-maximal mixing can also be produced by a dihedral flavor group with 
preserved subgroups. 
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1 Introduction 



Discrete groups have been widely used as flavor symmetry. However, only in some special cases 
there is a direct connection between the flavor group Gp and the resulting mixing pattern for 
the fermions, i.e. a correlation which does not rely on further parameter equalities not induced 
hy Gp- Such cases occur in the A/^ (T') models [2,3] as well as in our systematic study of the 
dihedral groups [1] where the key feature is the fact that Gi? is broken in a non-trivial way, i.e. 
one has to demand that certain of its subgroups arc preserved in different sectors of the theory. 
Especially, the fact that sizable mixing results from the mismatch of two different (directions of) 
subgroups is used in the A4 (T') models as well as in an application of the group D7 shown in [1]. 
In the group A4 {T') which has been studied in great detail tri-bi-maximal mixing (TBM) in the 
lepton sector is predicted, if one assumes that the left-handed leptons transform as a triplet under 
^4 (T'), and the left-handed conjugate leptons, e'^, and t^, as the three non-equivalent one- 
dimensional representations of the group. There exist two sets of gauge singlets which transform 
non-trivially under ^4 (T'): one set only couples to neutrinos at the leading order (LO), while the 
other one only to charged leptons (fermions). The flrst one breaks A^ (T') spontaneously down to 
Z2 {Zji) and the latter one down to Z3. The lepton mixing then stems from two sectors in which 
different subgroups of A4 (T') are conserved. In contrast to this, the up quark and down quark 
mass matrix preserve the same subgroup at LO. Similarly, it has recently been shown that such a 
mechanism can also be implemented with other discrete groups, for example the dihedral groups 
Dn and D[^. In a first application we observed in [1] that the Cabibbo angle 9c or equivalently the 
CKM matrix element | Vus \ can be predicted in terms of group theoretical indices only, such as the 
index n of the group D„, the index j of the representation under which the (left-handed) quarks 
transform and the misalignment of the two different (directions of) subgroups Z2 =< BA"*" > 



There is a crucial difference between the models using a dihedral symmetry and A4 (T') as flavor 
symmctr}^ namely the issue whether the representations under which the Higgs (flavon) fields 
transform are chosen or not. In our study on dihedral symmetries [1] we always assumed that 
for each representation fj, which (has a component which) transforms trivially under the relevant 
subgroup there exist(s) (a) Higgs field(s) transforming as fj, and acquiring a non-vanishing vacuum 
expectation value (VEV). Due to this the resulting mass matrices are only determined by the 
choice of the fermion representations, the dihedral group and the preserved subgroups, but not 
by the choice of the scalar fields. This makes our results less arbitrary. However, in the case of 
the A4 (T') it is necessary to choose the transformation properties of the scalar fields properly, 
i.e. one has to exclude scalars which transform as non-trivial singlets under A4 (T') and couple 
to neutrinos at LO, in order to arrive at the TBM scenario [3,4]. 

In this paper we investigate the idea of [1] by constructing a viable (low energy) model for the 
quark sector. The gauge group is chosen to be the one of the Standard Model (SM). We study the 
mass matrices numerically in order to demonstrate that all quark masses and mixing parameters 
can be accommodated. We discuss the Higgs potential under the assumption that all involved 
fields are copies of the SM Higgs doublet. Furthermore, instead of accommodating all quark 
mixing angles at LO it is also worth studying setups in which the Cabibbo angle is predicted 
in terms of group theoretical quantities, while the two other mixing angles are zero. This can 
be done in at least two ways: a.) one can choose the representations under which the scalar 
fields should transform or h.) one can look at cases in which the preserved subgroup in each 



and Z2 =< BPJ^'^ >: 




(1) 
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sector is not only a Z^-, but rather is itself a dihedral group -D^, q > \. Finally, we motivate 
possible extensions of the model to the lepton sector by performing a detailed numerical study. 
Additionally, we show that nearly TBM can be also accommodated by using a dihedral flavor 
symmetry. 

The paper is organized as follows: in Section [2] we review the findings of [1] which we explore 
in more detail; Section [3] treats the mixing matrix Vckm only - in an analytic way as well as 
numerically; in Section S] we study a model for the quark sector which incorporates the idea 
presented in [1] and show that it accommodates both quark mixings and masses; in Section [S] the 
Higgs potential, belonging to one of the models of Section [H is discussed and a numerical analysis 
proves that the advocated VEV structure can be achieved. Section E] is devoted to ansatze in 
which only the Cabibbo angle is generated at LO. In Section [7] we perform the same analysis, as 
for the quark mixing matrix Vckm in Section [3l for the lepton mixing matrix Vmns in order to 
see whether the fact that one element of the mixing matrix is described only in terms of group 
theoretical quantities is also applicable here. Thereby, we assume that the neutrinos are Dirac 
particles as all the other fermions and possess the same mass ordering, i.e. that they are normally 
ordered. Finally, we summarize our results in Section [8l Appendix El contains the possible forms 
of the mixing matrices Vckm and Vmns-, in Appendix [B] the group theory of D-j, i.e. the flavor 
group used in Section [3] and Section [5l is presented. Further details of the study of the Higgs 
sector are delegated to Appendix O 



2 Basics 



In this section we repeat the findings of [1] concerning the possible structure of (Dirac) mass 
matrices with a non- vanishing determinant. They are of the form: 



Ml = B , Ma = B (2) 
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Ms = B C , (3) 

\ Q D E I 

M4 = I C D E \ and M^=\ B D E (4) 

L>e-*'^j ^e-^<^j / \ 5e-*'^j Eq-''^'^^-^^ d e'' ) 

where A^ B, C, D, E are complex numbers which are products of Yukawa couplings and VEVs, 
(j) = ^ m [n: index of the dihedral group, m: index of the breaking direction) and j,k are 
indices of representations. Regarding M4 notice that we presented in [1] the transpose of this 
matrix. However, a transposition in general only corresponds to the exchange of the transfor- 
mation properties of the left-handed and left-handed conjugate fields under the flavor symmetry 
and therefore does not change the group theoretical part of the discussion about the preserved 
subgroups. Furthermore, these matrices are determined up to permutations of columns and rows 
which also only correspond to permutations among the three generations of the fields. As in [1] 
we work in the SM and with the assumption that all Higgs fields H in the model are copies of 
the SM one. Therefore the displayed mass matrices are those for down-type fermions, i.e. down 
quarks and charged leptons. The corresponding ones for up-type fermions, i.e. up quarks and 
(Dirac) neutrinos, require some changes due to the fact that only the conjugates of the Higgs 
fields, eH*, couple to up- type fermions and we use complex matrices for the two-dimensional 
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representations of Z?„. According to the rules of [1] on how to deduce the up- type fermion mass 
matrices from the shown ones, M4 and M5 are of the form 

Ce*'^'^ C 

M4=( Ce^'^j L>e*«^j E^e^'^j | and M5 = ( Be''^'^ De''t'^^+^'^ ^gi</.(j-k) | (^5) 

E D 

An exphcit example is given in Section [U where a model for quark masses and mixings is pre- 
sented. 

We concentrate on the last two forms, M4 and M5. This we do for two reasons: first we want to 
accommodate all masses and mixing parameters at tree level in the first part of the work, i.e. we 
do not want to rely on the fact that one mixing angle is only generated by higher order effects; 
second we would like to have the same mass matrix structure for up quarks (Dirac neutrinos) 
and down quarks (charged leptons). 

Let us briefly mention the origin of the matrix structures M4 and M5. The flavor symmetry is 
a single- valued dihedral group with arbitrary index n. The preserved subgroup is in both 
cases Z2 =< BA™ > where m = 0,l,...,n — 1. This subgroup allows non- vanishing VEVs for the 
following one-dimensional representations: 1^ (is always allowed to have a VEV), 1^3 for m even 
and I4 for m odd. All two-dimensional representations acquire a so-called structured VEV, i.e. 
for two fields ^/'i^2 transforming as an irreducible two-dimensional representation 2p their VEVs 

2 TT z p m 

have to have the correlation: < tpi >=< -02 > e . The notation of the representations 

used here is according to the one given in [1]. 

In case of M4 we take the left-handed fields L to transform as Ij^ ^^^^^ dihedral group, 
and the left-handed conjugate fields transform as the three singlets ~^—i2 is' ^ complete 
study of all possible assignments shows that one of the entries in the first row needs to be zero 
in order to prevent the determinant of the matrix from being zero. The matrix structure M5 
arises, if both left-handed and left-handed conjugate fermions transform as 1 -|- 2, L ~ (li,2j) 
and ~ (ll)2j^)- Here the constraint det(M) 7^ enforces the (11) element of the mass ma- 
trix to be non-zero, i.e. 1^ x l^j has to have a non-vanishing VEV. Note that the indices of the 
representations do not need to coincide, i.e. i 7^ 1 and j 7^ k is possible, although it might not be 
favorable from the viewpoint of a (partially) unified model. 

To study the mixing matrices arising from M4 and M5 for down-type as well as up-type fermions 
we observe that the products Mj m|, i = 4,5, can be written in the general form 

a be'l^ 
6e~*^ c 

where a, b, c, d and /3 are real functions of A, B, C, D and E. The phase [5 lies in the 
interval [0, 2 7r). Since we work in the basis in which the left-handed fields are on the left- 
hand side and the left-handed conjugate fields on the right-hand side, the unitary matrix which 
diagonalizes Mi mJ acts on the left-handed fields and therefore determines the physical mixing 
matrices, i.e. the CKM matrix and the MNS matrix. The three eigenvalues are given as {c — d), 
^ {a + c + d — y^{a — c — d)"^ + 8 6^) and ^ {a + c + d + y^(a — c — d)^ -|- 8 6^). Assuming this 
ordering of the eigenvalues the mixing matrix U which fulfills [/^ Mi mJ U = diag is of the form: 




U 



( cos(e)e^^ sin(e)e^^ 

]_J.<i>\ _sin(e) cos^e) 

y/2 \f2 

1 sin(e) cos(6>) -i^\ 

V2 V2 V2 
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The angle 9 is determined to be tan(2 0) = ■ Therefore it Hes in the interval [0, ^). If the 

three eigenvalues are not degenerate, the eigenvectors are determined by them up to phases 0- 
Therefore the variants of the mixing matrix U are given by permutations of the columns. With 
this at hand we can look for possible interesting structures in the mixing matrix which is just the 
product of two matrices of this form, as we assume that the up quark (Dirac neutrino) and the 
down quark (charged lepton) mass matrix is either of the form M4 or M5 . The mixing matrix is 
then of the form V = Wf VFJ with Wi being a variant of the matrix U above. For V = Vckm we 
have Wi = Uu which is the unitary matrix diagonalizing the up quark mass matrix and W2 = 
which is the corresponding matrix for the down quarks. In case oiV = VmnS: ^1 is equivalent 
to Ui and W2 to Ui, d. The matrix Wi contains the group theoretical phase according to the 
breaking direction rrii, the angle 6i and the phase For Wi = Uu we also use the notation (pu, 
niu, 0u and An analogous convention is used for U^, Ui and Ui,. It turns out that one of the 
elements is determined by the index j of the representation 2j under which two of the left-handed 
fields transform and the difference of the group theoretical phases 0i and 02 only. These phases 
do not depend on the values of the parameters A,B, ... (and therefore also not on a,b,c,d and 
f3), but only on the index n of the group Dn and the indices mi and m2 being the parameters 
that determine the subgroup to which the Higgs fields break D„ down. Therefore this element is 
determined by fundamental values of the model only and not by an arbitrarily tunable number. 
The actual form of (the absolute value of) the element is 

1 



1 _)_ g* i'f'i-'f'i)} 



cos((0i - (t>2) ^; 



vr 

cos(- (mi - 7712) j) (6) 



n 



Note that this value is only non-trivial, if mi 7^ m2, i.e. the (directions of the) subgroups which 
are preserved in the up quark (Dirac neutrino) sector and the down quark (charged lepton) sec- 
tor are not the same, i.e. only their mismatch leads to non-trivial mixing. This element can be 
traced back to the eigenvectors which correspond to the eigenvalue c — d m the up quark (Dirac 
neutrino) and the down quark (charged lepton) sector, i.e. the product of these two eigenvectors 
gives rise to the element ^ (l + e* j) . Therefore the ordering of the eigenvectors in the up 

quark (Dirac neutrino) and down quark (charged lepton) sector determines in which position of 
the mixing matrix the fixed element appears. As Vckm = U^, the fact whether the eigenvalue 
c — d \s associated with up, charm or top quark mass determines the row in which the element 
appears while the choice of the eigenvalue c — d to be either md, rUg or m;, determines the column. 
Analogously, the choice whether mg, m^ or m,- is given hy c — d determines the row, while the 
column in which the element appears is given by the fact which of the (Dirac) neutrino masses, 
mi, m2 or ms, is equal to c — d. 



In [1] we already mentioned that we can accommodate the CKM matrix element \Vus\ by 
cos(^) ~ 0.2225, i.e. by taking n = 7 and for example j = 3, m^j = 1 and m^ = 0. Here 
we show first which of the other elements of Vckm can also be accommodated well by the form 
|cos(^ {ruu - md)i)\. 

The elements of Vckm are precisely measured [5] 

n Q7QfiQ+0.00024 9979+0.0010 |"^Q«+0•09^ y ir,-3 

\Vckm\ = I 0.2271^°:°°}° 0.97296i°:S (42.21^0:10) x IQ-^ 

+0.32\ ^ in-3 (A^ ri+0.12^ ^ 1 n-3 n nnni nn+0.000034 



^.14+^:^1) X 10--^ (41.611^:^^) X 10--^ 0.999100 



'0.000004 



^ Since the eigenvectors should be normalized their length is fixed to one. 

^Throughout the paper we assume that the neutrinos are Dirac particles for simplicity. Therefore Vmns has 
the same structure as Vckm, i.e. there are no (additional) Majorana phases present in the lepton sector. 
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together with the Jarlskog invariant [6] Jcp = (3.08lQ }g) x 10~^. We restrict ourselves to values 
of n smaller than 30, since the index n of the dihedral group -D„ is correlated to its order and 
larger values of n correspond to larger groups. Enforcing n < 30 leads to a group order smaller 
than 60 which seems to be reasonable. Then we see that we can put the elements of the 1 — 2 
sub-block, i.e. \Vud\, |Ks|, and \Vcs\, into the form |cos(^ (m^ - m^) j)|. As \Vcd\ ~ |Ks| 
holds to good accuracy, also \Vcd\ can be described well by cos(A^). Furthermore \Vud\ ~ \ycs\ can 
be approximated well as cos(y^) ~ 0.9749 which points towards the flavor group -D14. Note that 
the value of \Vud\ as well as of \Vcs\ can be accommodated even a bit better with cos(^) « 0.9730. 
However, we do not use this, as it needs the group 1)27 which is a group of order 54 and therefore 
already quite large. Note that, even if \Vus\ is taken to be cos(^), there is no unique solution 
which flavor symmetry has to be used and to which subgroup it has to be broken, since for 
example taking j = 1, ruu = 3, rud = and n = 7 leads to |cos(^ (m^ — w^) j)| = | cos(^)| as 
well as j = 3, = 1, rud = and n = 7 and also j = 1, rriu = 6, md = and n = 14. As 
cos (^)l equals | cos(^)|, this allows us to deduce further possible values for j, m^, md and n 
like j = 1, niu = 0, md = 4 and n = 7. 

In the next section we study the cases \Vus\ and \Vcd\ equal to cos(^) and \Vud\ and \Vcs\ equal 
to cos(y^) in greater detail and thereby check whether we can always adjust the two other mixing 
angles and with the free angles 9u and 6d and also the Jarlskog invariant Jcp with the 
difference of the two phases (5u and fid- 



3 Analysis of Vckm only 
3.1 Remarks 

There are six possible forms for U which correspond to different identifications of the eigenvalues. 
However, the fact that <C rric ^ mj and <C rrig <C allows only three of them, as the 
eigenvalue \ {a + c + d— \J{a — c — d)"^ + 8 6^) is smaller than ^ (a + c + d + y^(a — c — d)^ + 8 6^) . 
Therefore, we will only vary the position of the eigenvector belonging to the eigenvalue c — d, 
while keeping the ordering of the two others fixed. The three different forms of the mixing matrix 
U are then: 



U 



u' 



u' 








sin((9)e*^ 




sin(0) 


cos(6») 








1 


sin(6») 


cos(e) 













sin(6l)e*^ 


sin{6») 


L p*</'j 


cos(6») 


V2 








1 


cos(e) 




V2 











sin{e) 


cos(e) 


L p«0j 


V2 








cos{e) -id,] 


1 


V2 ^ 




V2 



Combining them leads to nine distinct possibilities for the CKM matrix whose forms are dis- 
played in Appendix El Since we already mentioned that we want to concentrate on the 1 — 2 
sub-block we only need to consider the four possible combinations which involve the matrices U 
and U' as mixing matrices. 
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3.2 Numerical Study 



In this section we discuss the results of our fits to the CKM matrix where we assume that one 
of the matrix elements in the 1 — 2 sub-block is determined by group theory, as explained in the 
preceding section. There are three free parameters in the mixing matrix: 9^^^ and « = — Pd- 
We use these to fit the other two mixing angles Of^ and ^Is as well as the CP violation Jcp- 
The forms of Vmix presented in Appendix lAl show that two of the elements \ Vub\, \Vcb\, \Vtd\ and 
I Vis I are determined by cos{0u,d) in each of the four different cases. As these elements are small, 
the free angles 9u and 9d are restricted to be 6d,u ~ f • Therefore 9d^u is expanded around ^, 
0d,u = f ~ ^d,uy (-d,u > 0. The resulting four CKM matrices are (up to the first order in e^^^) 



/ cos(f^) cos(^) cos(4f)ed \ 

\VSkm\ « cos(^) cos(^) i|(i + ef')ed-2e^"6„| (7) 

Vcos(^)e„ i|(l + ef 0e«-2e*"ed| 1 / 

/ cos(f^) cos(^) cos(^)ed \ 

I^CKA/I « cos(i^) cos(^) i|(l + c¥«)e,-2e-ej (8) 

V i|(l + e^0e„-2e'"ed| cos(^)e„ 1 / 

/ cos(^) cos(^) 1 |(i + e^0ed-2e'"eu| \ 

\VSkm\ « cos(i^) cos(^) cos(^)e, (9) 

Vcos(^)e„ i|(l + e¥^')e„-2e'"ed| 1 / 

/ cos(^) cos(^) i|(H-ef ■Oed-2e^"e„| \ 

\V^^km\ « cos(4f) cos(^) cos(^)e, (10) 

V i|(l + ef ^)e„-2e'"ed| cos(2^)e„ 1 / 

Without loss of generality we have set the representation index j to 1, the group theoretical phase 
4>u to zero {rriu = 0) and the phase (f>d to ^ {nid = 1, n = 14) for Eq.([7j) and Eq. (fTOll . while we 
take it to be ^ (m^ = 3, n = 7) for Eq.® and Eq.Q. 

Comparing Eq.([7|) to the best fit values of \Vub\ and \Vtd\ given in [5] leads to ~ 0.0366 and 
€d ~ 0.0178. The phase a is then mainly determined by the values of \Vcb\ and \Vts\- A numerical 
computation leads to a best fit for a ~ 4.810 El- Furthermore one can calculate Jcp in this case: 

Jcp = ^ sin (^y) sin (^^) sin (2 9d) sin (2 9u) sin (j^ - « 

~ t: sin sin sin - a) e„ 



2 V7/ V14y V14 

A similar analysis can be carried out for the three other matrices Vckm^ ^ckm ™^ ^ckm '^ith 
similar results which we have collected in Tabled) The value of Jcp belonging to V^j^j^f, i.e. 
J^p, is of the same form as J^^p- For V^j^j^^ and V^]^j^j one finds 

Jcp = Jcp = sin ( ^ ) sin ( ) sin(26'd) sin(26'u) sin ( 



1 . / 67r\ . / 37r\ . /Svr 



As one can see in Tabled! eu,d have to be larger in case of V^j^j^^, since they are determined 
by |V"cfe| and \ Vts\- In this way the expansion of 9u,d around f gets worse and the second order 
in eud becomes important. This can be seen best in \Vus\ ~ 0.2225 and \Vcd\ ~ 0.2225 which 



^We performed a fit of Jcp and all elements of |VcifM| excluding the one which is fixed by group theory. 
Instead of taking the (very small) experimental errors we simply assumed 10% errors for all quantities. 
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Parameters 


^CKM 


^ CKM 


^ CKM 


^CKM 


e« 


+0.0364 


+0.0427 


+0.00831 


+0.188 


erf 


+0.0177 


+0.00405 


+0.0433 


+0.191 


a 


4.810 


2.355 


1.764 


0.2056 



Table 1: Fit results for eu,d {du.d) and the phase a for Vckm with either \Vud\, |K, 
\Vcd\ or \Vcs\ being group theoretically determined. 



are lowered to 0.2186(5) such that the discrepancy between the experimentally measured value 
and the result of the fit gets larger. However, corrections from higher-dimensional operators and 
explicit breakings of the residual subgroups can lead to further contributions allowing all data to 
be fitted successfully. 

4 Analysis of the Quark Sector 

After having shown that one element of Vckm can be explained in terms of group theoretical 
indices only and studying this issue numerically we want to go a step further and construct a 
viable model at least for the quark sector which includes this issue. The model is viable, if we find 
a numerical solution which accommodates not only the mixing parameters contained in Vckm, 
but also the quark masses. Due to the strong hierarchy among the quarks this is a non-trivial 
task, although the number of parameters in the mass matrices and exceeds the number of 
observables. Furthermore we have to show that a Higgs potential exists allowing us to realize the 
desired VEV structure. In the simplest case we assume that all Higgs fields are SU{2)i doublets 
as the Higgs field in the SM. 

4.1 D'! Assignments for Quarks 

Here we present ways to produce the two matrix structures M4 and M5 shown in Eq.([3j) and 
Eq.dS]) with the help of the dihedral group Dj. Choosing Dj as flavor symmetry leaves us the 
possibility of either determining \ Vus\ or \Vcd\ in terms of group theoretical quantities as cos(^). 

4.1.1 Matrix Structure M4 

For M4, we assign the quarks to 

Qi ~ ll , ( ) ~ 2i , u\,d'i^l2 , ni,3,43 ~ ll (11) 

under Dj. In this way we can generate the matrix structure found in Eq.([3]) and Eq.([5]) for the 
down as well as the up quarks. Thereby we assume that the theory contains Higgs doublet fields 
transforming as 1^ and 2^, which we call Hg and Hi^2- As the relation between the mixing 
parameters of Vckm and the group theoretical indices only arises, if the flavor symmetry Dj 
is broken down to a subgroup Z2 =< BA™" > by flelds which couple to up quarks, while it is 
broken down to Z2 =< BA™'* > with 7^ by fields coupling to down quarks, we need an 
extra symmetry to perform this separation. In the SM this can be achieved by a ^^""^^ symmetry: 

d'r ^ -df and ^ -Hf, Hf ^ -Hf (12) 
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while all other fields Qi, uf, and ^''^^ invariant under Z^^^\ Note that in principle also 
a Higgs field transforming as I2 under Dj could couple directly to the quarks. However, if this 
field acquires a non- vanishing VEV, its VEV breaks the residual Z2 generated by < BA™ >. 
Therefore Higgs fields ~ I2 either do not get a VEV or do not exist in the model at all. In 
both cases they are not relevant in the discussion of the fermion mass matrices. So, we deal with 
six Higgs fields coupling to the fermions, ~ (1]^,+1), ~ (2i,+l) and ~ {l-y,—!), 

^12 ~ (—1' ~1) under x Zg""^^ The matrices are of the form: 

/ yf (H^)* {H^Y \ f yf (Hf) vi (Hf) 

Mu = yf {Hfr yl mr Vl (HiY and Ma = yf {H^) yf {H^) yl 

\-y^{H'^r yimr ytmy I \-yi{Hf) yHHf) yHHf) 

where y"'*^ denote Yukawa couplings. The VEV structure is taken to be: 

{Hf'-)>0, {Hf) = {Hi)=Vd, {H^)=Vue-'-r and {H^)=Vue'-r 

with Vd > and Vu > 0. The VEVs are required to be real apart from the phase which is 
necessary for the correct breaking to the desired subgroup of Dj. 

The parameters A, B, ... shown in Eq.Q and Eq.@ can be written in terms of Yukawa couplings 
and VEVs: 

Au = yl (H^) , Bu = y2 (H^) , Cu = ytvue'''-^ , = yl Vue'^ , Eu = yl e"^ , 
Ad = yi{H'',), Bd = yi{H^), Cd = yivd, Dd = yivd, Ed = ylvd 

together with 0„ = (m„ = 3), = (m^ = 0) and j = 1, as the left-handed quark doublets 
of the second and third generation transform as 2^^. The preserved Z2 subgroups are generated 
by BA^ and B in the up and the down quark sector, respectively. As we have not fixed the 
ordering of the mass eigenvalues, the question which of the elements of Vckm is determined by 
group theoretical quantities to be cos(^) cannot be answered at this point. 



4.1.2 Matrix Structure M5 

For the case of M5, see Eq.Q and Eq.@, we can assign the quarks to: 

under Dj. We then need five Higgs fields for each sector, i.e. for the up and the down quarks. 
These transform as 

( ) ~(2i,+l), ( ) ~(22,+l) 
(11,-1), ( ) ~(2i,-l), ( ) ~ (22,-1) 

where we again assumed the existence of an extra ^^^^^^ symmetry. Under this z^'^"^) the quarks 
transform in the same way as in the example above, i.e. only the down quarks acquire a sign. 
The mass matrices are then in terms of Yukawa couplings and VEVs: 

/ yf (H^)* y| (Hf)" (H^)* \ ( yf (H^) y^ (H^) yf (Hf) \ 

Mu = yl {Hfy yl {h\Y yl {H^T and Md = yf {H^} yf {hi) yf {H^) 

\yUH^r vliH^Y vUh^r J \yUHf) yfiH^) yiiK) j 
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where y"'*^ denote Yukawa couphngs. The VEV structure is assumed to be: 
{Hf'-)>0, {Hf) = {H^)=v,, {hf) = {hi) = w,, 

(H^) =vue~^ , {H^) = Vue^ , {h\) =Wue~^ and {hl^) = w„e^ 

with Vd,u > and Wd^u > 0. As above we only consider real values for the VEVs apart from the 
phases which are required by the desire to break down to a certain subgroup of Di. 
Compared to the form of M5 given in Eq.Q and Eq.dS]) we see that the parameters A,B, ... are 
given by: 

Au = yl {H^) , Bu = yt Vu e"^ , C« = , D^ = y^ WuS'^ , Eu = yl {H^) , 

Ad = yf{H^), Bd = yivd, Cd = yivd, Dd = yiwd, Ed = yi{Hf) 

together with (pu = ^ {rriu = 3), 4>d = {rud = 0) and j = k = 1 for up as well as down quarks, 
since all generations transform as + 2-^ in this setup. Therefore the preserved subgroups in 
the up and down quark sector are again Z2 =< BA'^ > and Z2 =< B >. 

Note that the shown assignments are not unique, since it is also possible to use another two- 
dimensional representation instead of 2^ for the fermions. Obviously, then also the transforma- 
tion properties of the Higgs fields have to be changed accordingly. 

From the viewpoint of unification the second assignment in which the left-handed as well as the 
left-handed conjugate fields transform as 1 -|- 2 is more desirable. However in this case we need 
at least five Higgs fields transforming as 1^, 2j, 2j with i 7^ j in order to arrive at the matrix 
structure M5. As we have to separate the up quark from the down quark sector, i.e. have to 
have Higgs fields which either couple to up quarks or down quarks, we need at least ten such 
fields. Since we want to show the minimal model, we constrain ourselves to the first case, i.e. 
matrix structure M4, in the following numerical study and the study of the corresponding Higgs 
potential and only give a numerical solution for the second matrix structure M5. 

4.2 Numerical Analysis of Quark Masses and Mixing Angles 
4.2.1 Matrix Structure M4 

Coming to our numerical results we take all VEVs to have the same absolute value 61.5 GeV 
which equals the electroweak scale 174 GeV divided by \/8; because our complete model includes 
eight Higgs fields 0. The Yukawa couplings are taken to be 

= 1.07967 • (-2-17704) yU ^ 2.55955 • e^d-^ieoQ) yU ^ ^ 9543 . . (2.43366) ^ 
y2 = 3.89557 • IQ-^ • (-2-28452) yU ^ 7 47229 • 10"^ • e^d-^^es) ^ 

yf = 2.52251 • 10"^ . (3-00267) ^ yd ^ ^ g^^^^ . ^q-2 . (-2.29202) yd ^ g 20874 • 10"^ • e^(-0-540i4) ^ 
yi = 8.95471 • 10'^ • (-2-13972) yd ^ ^ ^49^7 . -^^-4 . (-1.59912) 

The values of the quark masses are then 

TUu = 0.0017 GeV , rric = 0.62 GeV , nit = 171 GeV , 
nid = 0.003 GeV , m, = 0.054 GeV , = 2.87 GeV 

■^The additional two Higgs fields which do not couple to the fermions directly, are necessary in order to break 
accidental symmetries present in the Higgs potential which we discuss in Section [5l The equality of the VEVs is 
motivated by our numerical study of the Higgs potential which clearly prefers solutions in which the VEVs are of 
the same order, otherwise severe fine-tunings of the parameters in the potential are necessary. However, this does 
not exclude in general the possibility that for example mt <^ nit could be explained by a hierarchy among the 
VEVs of the Higgs fields coupling only to up quarks and those coupling only to down ones. 
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which correspond to the values given at Mz [7]. For VckMi we find: 



/ 0.97492 0.2225 3.95 x lO^^ 

\Vckm\ = 0.2224 0.97404 42.23 x lO^^ 

\ 8.11 X IQ-^ AIM X 10^3 0.9991 

and Jcp = 3.09 x 10~^. All these values are within a 10% error range [5]. Furthermore \Vus\ is 
given by cos(^) = 0.2225. Due to the ordering of the eigenvalues the mass of the strange as well 
as the one of the up quark is determined by -v/2|Crf| and v^ICmI, respectively. They therefore 
correspond to the eigenvalue (c — d) in the language of Section [21 

The Yukawa couplings lie in the range 10~^...l due to the strong hierarchy of the quark masses. 
However this can be explained by the Froggatt-Nielsen mechanism [8]. The quarks transform in 
the following way: 

qpNiQi) = +1 , qFN{Q2,z) = +2 , qFN{dl^2;i) = qpNiuD = +i , qFN{u2^^) = -1 

under the additional U{1)fn symmetry. As usual we assume a gauge singlet with g_FAf(t^) = — 1 
which neither transforms under Dj nor under ^^""^^ and which acquires a VEV (t?) at a large 
energy scale El- According to the choice of the U{1)fn charges only the Yukawa couplings 
exist at tree level while all other couplings require some insertion of the field, i.e. become non- 
renormalizable involving some power of with A being the cutoff scale of the theory. One can 
then re-write the given Yukawa couplings y^''^ in terms of new couplings y"'*^ and e = = 3- 10"^ 
as yj'''^ = y^''^ where x is determined by the charges of the quark fields under the U{l)F]\f. The 
values of y^' are then all of natural size: 

y^ = 1.07967 • (-2-17704) ~u ^ 2.55955 • e^d-^ieoQ) ~u 

y2 = 1.29852 • (-2-28452) ~u ^ 2.49076 • e^(^-2469) , 

yf = 0.840837 • (3-00267) yd ^ ^^^^^ . (-2.29202) -d 

yi = 0.099497 • e^ (-2-13972) ~d ^ 0.116574 • e^ (-1-59912) 

4.2.2 Matrix Structure M5 

For the second matrix structure M5, we also performed a numerical study with the mass matrix 
structure given above and found for example the following possible values for the parameters 

■^u,d) -Bu,d: — 

Au = 40.40221 • 6^(0-185452) ^ ^ 0.238084 • (-2-99845) ^ 

Du = 0.420584 • (-3-13931) ^ Eu = 0.984542 • (-0-^49532) ^ 

Ad = 2.233447 • e^i-^-^^^^'^) , Bd = 0.051223 • (-3-05165) ^ (jd 

Dd = 0.058343 • e^(-2-4i4ii) ^ ^ 0.056221 • (-2-37708) _ 

All values are given in GeV. The phases (p^^d can be chosen to be (pu 
The quark masses are then: 

m„ = 0.0017 GeV , rric = 0.62 GeV , mt = 171 GeV , 
nid = 0.003 GeV , m, = 0.054 GeV , nib = 2.87 GeV 

^Here we assume that the (7(1) fat is broken explicitly in other parts of the Lagrangian so that no Goldstone 
boson arises when 1!} gets a VEV. 



= 0.723926 -e^ (2-43366) ^ 
= 0.68986 -e' (-0-54014) ^ 



= 117.4875 -e^ (-0-234118) ^ 
= 1.271448 -e* (-0-751605) ^ 

= ^ and (pd = 0. 
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and the absolute values of Vckm- 



/ 0.97489 0.2226 3.95 x lO^^ \ 

\Vckm\ = 0.2225 0.97401 42.23 x lO^^ 

\ 8.11 X 10-3 41.64 X 10-3 0.9991 / 

together with Jcp = 3.09 x 10^^. These values match the experimental results quite well. Note 
here that this time not \Vus\, but now \Vcd\ is given in terms of the group theoretical indices, 
i.e. \Vcd\ = cos(^) = 0.2225. S ince iV^tsleo^p ^ \^cd\exp^ \^cd\ is now 2% below its experimental 
value. This is due to the fact that the eigenvalue (c — d) introduced in Section [2] is given by 
rUc in the up quark and by nid in the down quark sector. These masses can be expressed in 
a simple way in terms of the parameters D^^ d and E^^d^ namely rric = — S^e-*'^'''^! and 
md = \Dd — Ede^'^'^^l with 0^ = <;/'d = and k = l|j as shown above. Also here the hierarchy 
among the parameters A^^d, B^^d, ■■■ which are products of Yukawa couplings and VEVs may 
not be explained by a hierarchy among the VEVs. Assuming that all VEVs have the same 
absolute value, i.e. GeV 55GeV|I|, an additional C/(l)i?Ar is responsible for the fermion 
mass hierarchy. The quarks have the U{1)fn charges: 

qFN{Q2,z) = +2 and qFN{dl) = +1 

and the other fields have zero charge. The parameter e = (??: U{\)fn breaking field with 
<1fn{'&) = — 1, A: cutoff scale) should be around 8 • IQ-^. 



5 Higgs Sector 

In this section, the Higgs sector belonging to the first numerical example given in Section [4.1.11 
is discussed. As already mentioned above, we concentrate on a multi-Higgs doublet potential. 
We are aware of the fact that such multi-Higgs doublet models usually suffer from the problem 
that large FCNCs are induced by the additional Higgs fields. However, as a proof of principle 
that we can produce our required VEV configuration the consideration of such a multi-Higgs 
doublet model seems to be reasonable. The minimal number of fields needed in order to produce 
the fermion mass matrices is 2 x 3, i.e. the model includes three Higgs SU{2)l doublets, called 
i/g and -fff 2' coupling to the down quarks and three doublets, and i^f 2' coupling to the up 
quarks. This separation is necessary, since the key point of the study lies in the fact that a sizeable 
mixing angle, like the Cabibbo angle, can only arise from preserved (non-trivial) subgroups, if 
we break to different (directions of) subgroups in the down and up quark sector. The subgroups 
correspond to different VEV configurations of the Higgs doublets {-fff,i?j^2} \^H^,Hf2\- 
An additional symmetry is introduced in order to perform the separation. According to 

Eq. p2|) in Section 14.1.11 the Higgs fields coupling to the down quarks acquire a sign under . 

We first construct the three Higgs doublet potential with Higgs fields Hg ~ 1^ and ( ^ ] ^ 



®k is here the same for the up and down quark sector due to the choice of the transformation properties of the 
left-handed conjugate fields M2,3 and 0(2,3. Note, however, that they could be in principle different. 

'^Here we assume the existence of only the ten Higgs fields which couple to the fermions in order to produce the 
mass matrix structure M5, as it seems unlikely that there are accidental symmetries in the Higgs potential which 
would enforce the existence of further Higgs fields. 
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The potential has the form: |f| 

1=1 \i=l I 



(14) 



A2 {h\ Rx - ul H2Y + A3 \hI H2\ 



+ ai {Hi H,) ij2HlHA+ {a2 {hI H^) {hI H2) + h.c.} + ^3 ^ 1^1 H,\^ 

\i=l / i=l 

As already shown in [9] and also mentioned in [10], this potential has an additional U{1) symme- 
try, i.e. there exists a further U{1) symmetry in the potential apart from the U{1)y symmetry. 
This further symmetry is necessarily broken by our desired VEV structure such that a massless 
Goldstone boson appears in the Higgs spectrum which is not eaten by a gauge boson. This 
problem cannot be solved by taking into account the whole potential for all six Higgs fields, since 
even if the terms, coupling the fields H^, i/f 2 ^12 together, are included, we find an 

accidental U{1) symmetry in the potential. Therefore we have to enlarge the Higgs sector by 
further Higgs fields in order to create new Dj invariant couplings which break this accidental 
symmetry explicitly. We find that this can be done in the simplest way by adding two Higgs 
fields transforming as 22 under D7. Due to their transformation properties they do not directly 
couple to the fermions (see Section I4.1.ip . We decided to add two such fields to the three Higgs 
fields which couple to the down quarks. Therefore the model contains eight Higgs doublet fields 
in total: three of them couple to up and three of them to down quarks, while the other two ones 
are needed for a viable Higgs sector: 

( ) ~(2i,+l), (15) 

(fj)~(2i,-l) and ( ) ~ (22,-1). 

under D-j x Zg""^^. The complete potential consists of three parts: 

V = Vu + Vd + Vmi^ed (16) 

where Vu denotes the part of the potential which only contains Higgs fields coupling to the up 
quarks, Vd contains the five Higgs fields which have a non-vanishing ^^""^^ charge (three of them 
give masses to the down quarks), while Vmixed consists of all other terms. The explicit form of 
the potential is given in Appendix O 

The VEV structure of the fields i/f and i/f '2 is determined by our desire to break down to two 
distinct Z2 subgroups in the up and the down quark sector (see Section I4.1.ip : 

{Hf'')>0, {Hf) = {Hi)=Vd, {H'^)=Vue-"-r and {H^)=Vue'-r 



with Vd > and Vu > 0. In contrast to this, the VEV structure of the fields X12 '^ot fixed in 
this way. However, in order to preserve the Z2 subgroup generated by B not only through the 
VEVs of the fields and Hf ^, but also by the VEVs of the fields xt2> ixf) = (X2) > ^ih be 
assumed in the following (see Section [2]). 

We proceed in the following way in order to find a minimum of this potential which allows for 



Note that (J2 is complex, but it can be made real by appropriate redefinition of the field Hs, for example. 
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our choice of VEVs: first we treat Vu and separately to find a viable solution for these two 
parts of the potential. Note that we can allow all parameters in the potential Vd to be real, as the 
VEVs of the corresponding Higgs fields are also real. Since Vu suffers from the above mentioned 
accidental U{1) symmetry, we find a fourth massless particle in the Higgs mass spectrum. In 
a second step we add as many terms as necessary from Vmixed to g^t a minimum of the whole 
potential V which does not have more than the usual three Goldstone bosons. It turns out that 
it is sufficient to take into account three terms in addition to F„ and Vd to get a viable solution. 
The terms are of the form: 



Note that we take all VEVs to have the same absolute value, since this considerably simplifies 
the search for a numerical solution, as a fine-tuning of the parameters in the Higgs potential is 
avoided. However, in principle other solutions should also be possible, e.g. the fact that the up 
quarks are much heavier than the down ones could be explained by assuming that the VEVs of 
the fields (much) larger than the ones of the fields i^f , -fff 2- 

Finally, let us mention that the resulting Higgs masses are usually in between 50 and 500 GeV. 
These values are either not favored by the constraints coming from FCNCs or already excluded by 
direct searches. There are two reasons for the too low Higgs masses: on the one hand Vu contains 
an accidental symmetry and on the other hand all mass parameters of the potential are chosen 
to be of natural order, i.e. to be around the electroweak scale. Additionally, all quartic couplings 
of the potential must be perturbative. However, as already mentioned above, this model is not 
intended to be fully realistic. Adding breaking soft masses to the potential might allow to 
push the masses of the additional Higgs particles above 10 TeV. 

The rest of the discussion of the potential is delegated to Appendix O where we present a nu- 
merical solution for the parameters of the Higgs potential and the resulting Higgs masses. 



6 Ways to generate 6c only 

In the preceding sections we confined ourselves to cases in which all mixing angles can be repro- 
duced at tree level. Therefore we only discussed the matrix structures M4 and M5 of Eq.dH) and 
Eq.([5]). However, and ^23 ^'^^ roughly an order of magnitude smaller than the Cabibbo angle 
9c = ^12 which gives reason for also considering matrix structures which lead to only ^c" 7^ at 
LO. For this a block matrix structure (with correlated elements), which we introduced in Eq.dSj), 
is suitable. Such a structure can be achieved in at least two different ways: a.) we can simply 
omit some of the Higgs fields which are in principle allowed a VEV in order to arrive at the 
zero elements of the mass matrix; h.) we can demand that the preserved subgroup is not just a 
Z2 symmetry, but a dihedral group Dq with g > 1 Note that due to the choice of the scalar 
fields in case a.) the results of such a model are a bit arbitrary, since the structure of the mass 
matrices and therefore the mixing pattern is not fully determined by the fermions and the flavor 
symmetry alone. In the following, we show examples for the two cases. For case a.) the simplest 
example is probably the one in which the quarks transform as 



The general results of the mass matrices from the various preserved subgroups can be found in [1]. 
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under D-j and we assume that there exist two sets of three Higgs fields transforming as 

(J;!)~(22,+1), ( J^^ ) ~(22,-l) 

under D-j x Zg""^^ one of them couphng to up and the other one couphng to down quarks. 
The additional z'^^^^ symmetry is the same as used above (see Section [4.1.ip . Then the mass 
matrices are of the form: 

Mu=\ ylWr Vsih^r and Md= {H^) yl {hf) 

V yf (if")- / \ yf (Hf) 

Assuming the VEV structure: 

{H^'n>0, {hi) = {hi)=wd, {h't)=Wue-"-r and (/i^) = «;„e^ 

with Wd^u > we arrive at 

/ y^wue"-r y^{H^) \ / yfwd yi{Hf) 

Mu= y^{H^) y^Wue^ and = yi{Hf) y^Wd 

V yf(F«) / V yf{H^) 

and 

/ |cos(f^)| |cos(^)| 0\ / 0.97493 0.2225 \ 

\Vckm\=\ |cos(^)| |cos(f^)| « 0.2225 0.97493 (17) 

\0 l/\0 01/ 

The preserved subgroup is a Z2 in each sector which is generated by B A'^ and B in the up quark 
and the down quark sector, respectively. The masses of the quarks are (m^, m^,m^) = (|?/2 + 
y^Wu\My^{H^)-y^wX,\yf (^^1') and (mlmlml) = {\yUHf) -yiwd\MyUHf) + yiwd\\ 
lyf (Hg)]"^), i.e. the mass of the third generation is solely determined by y^''^(i?"''^). Note that if 
the VEVs of Hs''^ are taken to be large in order to explain the large mass of the third generation, 
the Yukawa couplings have to be suppressed. This might be viewed as fine-tuning. A possible 
solution is the assumption of an additional U{1)fn as already used above or to consider the case 
b.) instead. The possibility to choose the two-dimensional representations for the left-handed and 
left-handed conjugate fields to be distinct from each other, such that the trivial representation 
1^ cannot be coupled to the first and second generation, does not exist, since in this case also 
two of the four zeros disappear. The reason can be found by looking at the Kronecker products 
shown in Appendix [BJ Actually, this setup is very similar to the one shown in Section 14.1.21 
The main difference is the fact that now there are no Higgs fields transforming as 2^ under Dj. 
The other difference is that the first and second generation of the left-handed and left-handed 
conjugate fields are unified into the doublet under the fiavor group instead of the second and 
third one as done above. However, this only leads to a change in the appearance of the mass 
matrix itself, but does not have any phenomenological consequences, since this permutation of 
fields is cancelled in the mixing matrix. The existence of six instead of ten Higgs fields coupling 
to the fermions may be advantageous with regard to the problem of FCNCs mediated by these 
fields. The corresponding Higgs potential ought to be of the same form as the one discussed in 
Section [5l 
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The second case b.) cannot be maintained with the flavor group -D7 which we used throughout this 
work, since it only contains Zq groups as subgroups, but no dihedral ones Dg, q > 1. Therefore 
we have to consider the group D14 instead. In the study of the Vckm elements in Section [2] and 
Section [3] Z)i4 turned out to be the smallest group which is appropriate to describe the elements 
\Vud\ and \ Vcs\ in terms of group theoretical indices. As argued in Section [2] and Section [3] it can 
also be used in order to reproduce the -D7 results, i.e. either \Vus\ = \ cos(^)| or \Vcd\ = \ cos(A^)|. 
Here we just show a possible example in which D14 is broken to its subgroup D2 =< A^, B A"^ > 
(m = 0, 1, 6) in order to reproduce a matrix of block structure. We assign the quarks to 

under D14. According to the Kronecker products 

1.1^^1 = 2i and 2]^ x 2]^ = 1 + I2 + 22 

the Higgs fields which can in principle couple to form L>i4-invariants have to transform as 1\, 
—2^ —1 —2- However, I2 is not allowed a VEV and the representation index j of 2j has to 
be even0. Therefore we take 

H^-ll, (fi)~22, Hf^li and ( f ^ ) ~ 23 

(with implicit ^^""^^ assignment as above) and arrive at the matrix forms which are exactly the 
same as given above for case a.) 0, if we assume the VEVs to be 

{K^'')>0, {H^) = vue-'r, {H^)=vue'-r, ^Hf) = {H^) = 

The subgroups which are preserved by the VEVs in the up and down quark sector are then of 
the form D2 =< A'^, B A™' > with m„ = 6 for the up quarks and rrid = for down quarks. Also 
here the Higgs fields Hs'^ couple to all three generations. In order to avoid this one can assign 
the quarks to different D14 representations, e.g. 

Since 2^ x 23 decomposes into 23 and 2^ in D14, the 1 — 2 sub-block of the mass matrices is 
produced by the VEVs of Higgs fields belonging to -D14 doublets instead of the sing let Hs' . As 
the indices of the representations 23 and 2^ are even, they are allowed a VEV by the requirement 
to preserve a D2 subgroup of -D14. We need five Higgs fields transforming as li + 22 + 24 for 
the down as well as the up quarks. The general form of the mass matrices reads 

/ VsiK)* Vi^H^Y \ / yiihi) yi{Hf) \ 

Mu = 2/2" {HfY and Md = {H^ yl (hf) 

V yl{H-Y I \ 

With the VEVs 

(i/f")>0, {Hf) = {Hi)=Vd, {hi) = {hi)=wd, 

{H^) =vue~^ , {H^) = u„e^ , (/i^) =Wue~^ and {h^) = Wue'~^ 

^"in general this index must be divisible by the group index of the dihedral subgroup which should be preserved. 
'^^The Clebsch Gordan coefficients necessary for the calculation of the mass matrices in D14 can be found in a 
general form in [1]. However, in this special case they coincide with those given for the group D7. 
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for Hs''^ ~ 1\, H^2 ~ —2 ^12 ~ —4' ^™ achieve 

y^Vue^ y^Wue-^ and = y^Vd yiwd 

yY {Hf) J V yf (H^) J 

For {ml,ml,m'f) = + «^«|My>« - ■w«|Myi and (m^.m^.mg) = {\y^Vd- 

yfi^dP) I2/2 + 2/3 f^dP) lyf (-f^s)P) the CKM matrix is of the form as given in Eq. (fT7|) . Although 
the Higgs fields Hg''^ couple in this setup only to the third generation and therefore can have 
a large VEV without spoiling the masses of the lighter quarks, there still exists a source of 
fine-tuning, since the uncorrelated parameters 2/2 3) ''^d,u and Wd,u have to be arranged such that 

1^2 "Vd - yi Wd\ <C 1^2 Vd + yi Wd\ for rud <. rUg and Vu + y^Wul <^ |y| Vu - ^3 for ruu <C mc- 
The preserved subgroups in the up and the down quark sector are again D2 =< A'',BA^ > and 
D2 =< A^,B >, respectively. 



7 Numerical Analysis of Vmns 

A similar analysis as done in the case of Vckm can also be carried out for the lepton mixing 
matrix Vmns- We assume that the neutrinos are Dirac particles as all the other fermions and that 
they have the same ordering as the other fermions, i.e. the neutrino mass spectrum is normally 
ordered. This allows us to use the matrix structures found in Appendix |A] also for Vmns- Since 
the entries of Vmns are not strongly restricted by experiments [11] (at 3 a): 

/ 0.79 -0.88 0.47- 0.61 < 0.20 \ 
I^mS I = 0.19 -0.52 0.42 - 0.73 0.58 -0.82 (18) 
V 0.20 -0.53 0.44 - 0.74 0.56 -0.81 / 

there are several more possibilities to accommodate the various matrix elements regarding the 
choice of the group index n, and the values mi, rrii, and j. However, as we intend to build a 
model which includes quarks as well as leptons, we stick to the selected values of n, n = 7, 
n = 14, which fit the CKM matrix elements of the 1 — 2 sub-block best, if we restrict ourselves 
to small n. We check element by element of Vmns whether we can put it into the form | cos(^)| 
where Z = 0, 1, 2, 6 or | cos(^)| with / = 0, 1, 2, 13. According to Eq.((TS]) all elements of the 
second and third row can be approximated by a cosine of the form | cos(4^)| (/ = 0, 1, 2, 6) or 
I cos(^)| {I = 0, 1, 2, 13) El. We take into account all possibilities shown in Table[2]and perform 
a numerical fit of the mixing angles ^12, ^13 and ^23- In the fit procedure we compute the sines 
of the three mixing angles and compare these to the best fit values, which are sin^(^23) = 0-^) 
sm [e\{) = 0.3 and sin^(0-i^2) — [12] r^ . Again, we replace the experimentally allowed 2cj or 3(T 
ranges by 10% ranges (around the best fit value). For sin^(0i3) we consider two possible upper 
bounds: sin^(0i3) < 0.025 which corresponds to the 2a bound [12] and a much more loose bound 
sin^(0i3) < 0.1 being even larger than the 4cr bound [12]. This is done, since the numerical study 
showed that loosening the bound on sin^(^i3) leads to several more solutions. Our results for 
sin^(^i3) < 0.1 are summarized in Table [3] where we also display the numerical values for 9i, O^, 
and a = Pi — Pu together with the resulting mixing angles and the (Dirac) CP phase 6. 

^^We omit the trivial possibility that the (13) element can be approximated by 0. 

'^^Note that these best fit values are not presented in the same global analysis as the above mentioned allowed 
3(T ranges for the elements of Vmns- Nevertheless the deviations are very small such that we do not consider this 
to lead to a major difference in our numerical analysis. 
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Element (ij) 




(21) 


cos(^ 


(22) 


cos(ff 


(23) 


cos(^ 


(31) 


cos(^ 


(32) 




(33) 


cos(^ 



Possible cosines 



)(^ 



-)(^ 
-)(^ 



0.2225), 
0.4339), 
0.6235), 
0.2225), 



cos 



cos( 
cos( 



COS I 
COS I 



3^ 
14 



-) 0.6235), cos 



:if)(^ 

0.6235) 

14 ^ V 



0.4339) 
0.6235) 
0.7818) 
0.4339) 

0.7818) 



Table 2: Possibilities for the group theoretically determined element in Vmns- 
Note that, e.g. cos(^) equals cos(S^), i.e. it could also be reproduced in the 
group Z)i4 with j = 1 and mi — rriv — 6 and not only in D7 with j = 1 and 
mi — mi, = 3. Furthermore, for example, cos(4^) is also included implicitly in the 
hst, as |cos(if )| = |cos(2^)|. 



One can observe the following: there are some cosines listed in Table [2] for which no fit with 
< 1 has been found. In all these cases the value of the fixed Vmns element lies almost outside 
the ranges shown in Eq. ()18p . e.g. for the (23) element the possible cosines are cos(^) « 0.6235 
and cos(^) ^ 0.7818 with the first being quite close to the lower bound (0.58) and the second 
one close to the upper one (0.82) of the allowed range. Furthermore, by having a closer look at 
the form of jV^/^^l given in Appendix El one realizes that tan(^23) is simply determined by the 
expression: 



tan(6'; 



23 J 



cot 



Mi 



COS (6// J 



cot 



005(6*/) 



(19) 



Taking the argument of the cotangent to be either ^ or ^ leads to the numerical values 



tan(6'; 



23) 



0.7975 



005(6*;) 



or tan(6l23) ^ 1.254 



cos{9i) 



At the same time the sine of 9i is determined by the 

I - M i 

0.7818 I sin(6'z)| and 0.6235 | sin 



r23 



of sin(6'i3): liVj^f^gjis 



^and^ |(y|3 



sm 



sm 



13) element of Vmns, i-e. by the value 



sm 



TT {mi-m„)i 



which gives for 

MNS)'^3\ = u. (Olo I sm(^p;j| anu u.uzoo | 5un^(7;j|, respectively, i.e. | 5in(0/)| has to be 
as small as possible to fulfill the experimental bound on sin^(^i3). Then | cos(0;)| ~ 1 holds so 
that we can deduce the approximate values 0.7975 and 1.254 for tan(023) from Eq.([T9]). These 
correspond to 6*23 ~ 38.57° (sin2(6'23) ^ 0.3888) and 023 ~ 51.43° {sin'^{923) ^ 0.6113), i.e. they 
are at the boundaries of the 2 a range for sin^(023) [12]. Similar statements hold in case of |V^/7vsl- 
Furthermore, one observes that in all cases the CP phase 6 is trivial, i.e. or vr with a numerical 
precision of 0(10"^). Therefore Jcp always vanishes. In order to understand this result, we have 
a look at the formulae given for V^^^, V^^L' ^mix ^mix Appendix |Al As a common feature 
the (13) element of the mixing matrix is given by 



1 



[-(1 + e- 



-i ict>i-(t>i,)} 



) sm{ei) cos{e^) + 2e^" cos{9i) sm{9^)] 



(20) 



In all cases, 9i and 9i, are predominantly determined by one element of the first row and the third 
column of Vmns, respectively. Then a can be used in order to minimize the absolute value of 
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Element 


Cosine 


ei 


Ou 


a 


sm2(0i2) 


sin2(023) 


sin2(0i3) 


5 


(21) 


cos(^) 


0.9790 


0.7881 


4.937 


0.2957 


0.5085 


7.037 X 10-2 


~ vr 


cos(^) 


1.1829 


0.6725 


5.161 


0.3001 


0.4999 


6.173 X 10-3 


~ 


(22) 


cos(^) 


- 


- 


- 


- 


- 


- 


- 


cos(^) 


0.7728 


0.4486 


5.386 


0.2999 


0.4996 


6.668 X 10-3 




(23) 


cos(^) 


- 


- 


- 


- 


- 


- 


- 


cos(4f) 


- 


- 


- 


- 


- 


- 


- 


(31) 


cos(Af ) 


0.9790 


0.7881 


4.937 


0.2957 


0.4915 


7.037 X 10-2 


~ 


cos(^) 


1.1829 


0.6725 


5.161 


0.3001 


0.5001 


6.173 X 10-3 


~ TT 


(32) 


cos(^) 


0.7728 


0.4486 


5.386 


0.2999 


0.5004 


6.668 X 10-3 


~ 


(33) 


cos(^) 


- 


- 


- 


- 


- 


- 


- 


cos(^) 

















Table 3: Numerical results for Vmns in case of sin2(6»i3 ) < 0.1 and 10% errors for 
the other two sine squares. 5 is given with a precision of OilQ^''). 



the (13) element of Vmns- A minimization with respect to a shows 

j TT 

a = -{(t)i - (t)v) - + TT y = {mi-m^)i + 'Ky with y € Zq (21) 

2 n 

The minimum value for I sin(0i3)| is then I cos(((/);— (/),y) ^) sin{9{) cos(0,y)+(— l)^"*"-*^ cos(6';) sm.{9u)\- 
However, in all cases the expression is only minimized for y = 0, 2, as the involved sines and 
cosines are all positive (remember that 6i and Oy are restricted to be smaller than ^ by definition 
and also {(pi — 2 — n ~ '^i^) j which is the argument of the cosine displayed in the tables is 
always smaller than ^). As Jcp is proportional to sm{{(j)i — (pu) 2 +'^)' it is zero for the calculated 
value of a. Therefore 5 must be either or vr. Additionally, we found an explanation for the val- 
ues of a shown in TableOgiven in terms of the group theoretical quantities, i.e. 2 vr — ^ ~ 4.937, 
2 7r — -^j ~ 5.161 and 2 7r — ^ ~ 5.386. Since a has to lie in [0, 2 7r), y equals two in all cases, 
see Eq.dnj. 

As a last observation we report that there exist similarities among the different cases, e.g. fixing 
the (21) element to be cos(^) is similar to fixing the (31) element to the same value. In both 
cases the fit values of 61, 9y and a are the same. Therefore, the results for siv?{9i2) and siv?{9i'^) 
coincide (up to O(10~^)), while sin2(023) is shifted from being 0.5 + e to 0.5 — e with e « 0.0085 
and the CP phase 5 shifts from vr to 0. Looking at the mixing matrices one recognizes that these 
similarities are due to the fact that the second and the third row are interchanged. 
Using the 2a bound s\v?{9i2,) < 0.025 no solution with < 1 is found in the cases in which 
the (21) or the (31) element is fixed to the value cos(^), since the values for s\v?{9i'i) shown in 
Table [3] are quite large. For the other configurations we again find viable fits in which the values 
9i , and a are very similar to the ones given in Table El 

Until now we only investigated the cases in which the group theoretically fixed element is given 
by one of the cosines shown in Table [21 However, as already remarked several times we can also 
look at cases in which the cosine is for example cos(^) instead of cos(^), since cos(^) is just 
the negative of cos(^). In terms of group theoretical quantities this corresponds to sending mi 
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Element (ij) 



Possible cosines 



(11) 

(12) 
(21) 
(22) 
(31) 
(32) 



cos 

COS 



COS I 
COS I 



cos( 
cos( 



)(^ 



0.7818) 
0.6235) 
0.4339) 
0.6235) 
0.4339) 
0.6235) 



Table 4: Possibilities for the group theoretically determined element in Vmns, if TBM 
is assumed to be the best fit. For further conventions, see Table [2] 



to n — mi and therefore (/>/ to 2tt — (j)i 0. The general forms of the mixing matrices given in 
Appendix \K\ show that such a transformation does not change the absolute values of the matrix 
elements, if we replace the phase a by —a (2 tt — a) at the same time. In contrast to this Jcp is 
not invariant and changes its sign. In the analysis of the leptonic mixing parameters this is not 
relevant, since the phase(s) have not been measured. Moreover, in all cases considered here Jcp 
is almost zero (up to 0(10"^)). Therefore, we get the same results for these equivalent cases. 
Note, that in case of the quark mixing matrix we would have to expect different results, since 
there Jcp is known from experiment and its sign change leads to a distinct solution. 
Apart from studying how well one can accommodate the experimentally allowed ranges, it is also 
interesting to see whether one can reproduce some special mixing pattern in the lepton sector. 
In the following we discuss the TBM scenario which has initially been discussed in [13], since all 
elements of the lepton mixing matrix can be written in terms of fractions of square roots 

and^: 



^tTBM 
^MNS 



2 1 

^/^ 




(22) 



corresponding to sines of the mixing angles: 

sin2(0™) = i, sin2(0™) = i and sin2(e™) = . 

However, it turned out to be not just an assumption of a special form of Vmns-, but it is a 
robust outcome of certain models based on the discrete non-abelian symmetries A/^ or T' [2,3]. 
Therefore we want to analyze whether we can also accommodate the TBM with mixing matrices 
of the form Vmix as given in Appendix|Xl The uncertainty in the mixing matrix elements is taken 
to be 10%, i.e. the fixed element given by cosine | cos(^)| for I = 0, 1,2, ...,6 or | cos(^)| with 
/ = 0, 1, 2, 13 should lie in one of the ranges: 

/ 0.73 - 0.90 0.52 - 0.64 < 0.20 \ 
^mS^ = 0.37-0.45 0.52 -0.64 0.64- 0.78 (23) 

V 0.37 - 0.45 0.52 - 0.64 0.64 - 0.78 / 

The bound on the (13) element is taken to be the same as in Eq. (ll8p . As shown in Table [H the 
elements (11) and (12) can now be described by a cosine of the announced form, while we find 



'^^Thereby, we have set <pf to zero and j to 1 without loss of generality. 
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Element 


Cosine 


Oi 


eu 


a 


sin2(0i2) 


sin2(023) 


sm2(%) 


5 


(11) 


cos(^) 


0.4396 - 1.131 


1.139 


e [0,27r) 


0.3441 


0.5000 


6.808 X 10-2 


G [~ 0, ~ 27r) 


(12) 


cos(^) 


- 


- 


- 


- 


- 


- 


- 


(21) 


cos(ff) 


1.132 


0.6697 


5.161 


0.3334 


0.5000 


1.968 X 10-3 


~ 


(22) 


cos(^) 


0.8235 


0.4557 


5.386 


0.3331 


0.4991 


1.245 X 10-2 


~ vr 


(31) 


cos(^) 


1.132 


0.6697 


5.161 


0.3334 


0.5000 


1.968 X 10-3 


vr 


(32) 


cos(^) 


0.8235 


0.4557 


5.386 


0.3331 


0.5009 


1.245 X 10-2 


~ 



Table 5: Numerical results in the case of TBM. We assume that the bound on sin'^ (^13) 
is 0.1 and 10% errors for the other two sine squares. The values of 5 have a numerical 
precision of 0(10"''). Note that in case of the (11) element being cos{^) 5 can take 
arbitrary values, (for details see text). 



less possibilities for the other elements compared to the case of the experimentally allowed range, 
see Table El This analysis is analogous to the one above. Again, we display the results for the 
fits using the loose bound for s\v?{6i'i), sin2(6'i3) < 0.1. Similar to above, there is a case in which 
we have not found a fit with < 1. For the cases in which either the (21), (22), (31) or (32) 
element is determined by group theory all statements made above can also be applied here, i.e. 
the CP phase 5 is either or tt, the phase a is fixed to a certain value which minimizes | sin(0i3)| 
and there exists a similarity among the cases with a fixed (21) ((22)) element and a fixed (31) 
((32)) element. Therefore, we focus on the discussion of a group theoretically determined (11) 
element of Vmns- This case exhibits some new features not present in the other ones. First of 
all, we find that 9i can take values in a certain range instead of being fixed to a single value. 
All of them lead to the same mixing angles. The same is true for a which varies between and 
2 7r. This is related to the fact that we do not fit the CP phase 5 (or equivalently the Jarlskog 
invariant Jcp)- As a result Jcp can take any value in the range (-5.776.. .5.776) x 10-^. We 
observe that 9y is fixed by the fit of si'c?{9i2) and siT?{9\'j,). Fitting them at the same time leads, 
unfortunately, to a too large value for s\v?{9i2,) (see Table [5]). The allowed range for 9i can then 
be found analytically under the assumption that sin2(023) = \^ since in this case the (23) and 
(33) element of Vmns have to be equal. Equating the expressions |(V'^L)23p and |(l^mL)33p 
found in Appendix |X] leads to 



tan(2.,)= s.n^(^.)-co.^((0,-0.)i)co.^(^.) 

C0S(((/); -(t)u)^ + a) cos((0z - (pu) ^) SUl{29u) 

with 9i, determined by s\v?{9i2.i2,)- Allowing a G [0, 2 7r) one finds the maximal range of 9i to be 
z < 9i < ^ — z with z ~ 0.4396 for 9y ~ 1.139 and {(pi — 4>,y) 2 = ij which corresponds to 
the numerical values given in Table [5l Furthermore, Eq. p4p shows that 9i is a function of a. 
Demanding sm'^{9i3) < 0.025 removes the possibility that the (11) element of Vmns is determined 
by group theory, while it leads to expected slight changes in the results of the fits for the rest of 
the cases. As expected, a comparison of these results to the ones for the experimental best fit 
values shows that there are only small changes in the precise values of 9i, 9y and a. 

In this section we have shown that it is possible to fit the lepton mixing angles [12] in a framework 
in which one of the elements of Vmns is completely determined by group theoretical quantities of 
a dihedral fiavor symmetry. We restricted ourselves to the dihedral groups D^ and -D14, since they 
allow us to explain the Cabibbo angle via group theory. A main result of the analysis is that Jcp 
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vanishes in all cases. Using the formulae given in Appendix |X] one can show that the vanishing 
of Jcp is correlated with the minimization of |sin(0i3)|. As the bound on |sin(0i3)| is the 
strongest constraint on the solution, we showed our results for a very loose bound. Furthermore, 
we analyzed how well one can mimic the TBM scenario. We found several possible solutions. 
One of these is of special interest, since it also allows for non-trivial CP violation. However, 
the corresponding value of sin(0i3) is very large and therefore this solution is disfavored. These 
results demonstrate that it is possible to treat the lepton mixings in the same way as the ones of 
the quarks. Small corrections are expected in a complete model, e.g. due to explicit breakings 
of the preserved subgroups. 



8 Summary and Conclusions 

In [1] we studied the dihedral groups as possible flavor symmetries. The key feature there is the 
fact that the flavor symmetry is not broken in an arbitrary way, but one requires that a subgroup 
has to be preserved in all cases. It turned out that the number of possible mass matrix structures 
which arise from such a setup is very limited, if we assume that the mass matrix has a non- 
vanishing determinant. As a first application we discussed in [1] the possibility to describe one 
element of the CKM mixing matrix only in terms of group theoretical quantities, i.e. the index n 
of the dihedral group Dn , the index j of the representation 2j under which the fermions transform 
and the indices and nid of the residual subgroups Z2 =< BA™"" > and Z2 =< BA™*^ >: 



1 _l_ g* {<l>^-'t>d)i 



COs{{(t>u - 4>d) ^) 



cos(- (m„ - md) j) (25) 
n 



where = ^ mu and (f>d = ^ n-d- Eq. ()25l) shows that a non-trivial mixing angle demands 
iTT-u / rrid, i.e. the two Z2 subgroups have to be distinct. It has been pointed out that |T4s| can 
be fitted well with cos (^) 0.2225. In this work, we first studied which of the other elements 
of VcKM can also be approximated well by Eq.(|25p for certain values of the group index n. For 
the smallest two appropriate values of n, n = 7 and n = 14, each element of the 1 — 2 sub-block 
of VcKM can be put into this form, i.e. I^^^l « \Vcs\ ~ cos(yj) and \Vus\ ~ iKdl ~ cos(^). A 
numerical analysis showed that the other two mixing angles, Of^ and 02^, and the CP phase 6 can 
be fitted well with the free angles and the phase a = f3u — Pd- Since the fixed element cannot 
be fitted, the results for Vckm are very close to the experimental values, but not within the 
(very small) experimental errors [5]. However, several sources of corrections exist in a complete 
model, e.g. possible higher-dimensional operators as well as small, but explicit, breakings of 
the preserved subgroups. In a next step, we presented a low energy model for the quark sector 
which incorporates the described idea. The flavor symmetry is taken to be D^. It is broken only 
spontaneously at the electroweak scale by Higgs fields transforming as doublets under SU{2)l. 
With a numerical fit we showed that all quark masses and mixing parameters can be fitted well at 
the same time. As the VEV configuration determines the subgroup to which the flavor symmetry 
is broken, it is necessary to investigate whether this can be achieved by the Higgs potential. We 
studied this issue for the minimal model in which the three Higgs fields Hg and H]^2 couple 
to up quarks, while the fields Hg and iff 2 couple to down quarks only. Unfortunately, the 

Higgs potential containing only Hs'^ and Hf'2 has an accidental symmetry which is necessarily 
broken by the desired VEV configuration. For this reason we had to add two further Higgs 
fields to the model. These do not directly couple to the fermions due to their D7 transformation 
properties. The couplings of Hs'^ and iff '2 to these Higgs fields break all accidental symmetries 
of the potential. A numerical study showed that the needed VEV configuration can be achieved 
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with this potential. However, there are two obstacles: first of all if the quartic couplings of the 
Higgs potential are in the perturbative range and the mass parameters are taken to be around 
the electroweak scale, the Higgs masses turn out to be too small, i.e. some of them are even 
below the LEP bound [14]. This could be cured by adding mass terms which break the flavor 
symmetry softly in the Higgs potential and allow for larger Higgs masses. However, even then 
this model might suffer from the problem that FCNCs induced by the additional Higgs fields are 
too large to pass the experimental bounds. The second obstacle is the fact that we are only able 
to accommodate the VEV configuration as one possible solution of the Higgs potential, but not as 
a favored solution. Moreover, there is in general no way to stabilize such a configuration against 
further corrections in a multi-Higgs doublet potential. Therefore this model is meant as a proof 
of principle rather than a realistic model. A way to circumvent these problems is to disentangle 
the scales of the electroweak and the flavor symmetry breaking by using flavored gauge singlets 
instead of Higgs doublets and thereby break the dihedral symmetry at higher energies [4] . 
Accounting for the fact that the Cabibbo angle 9c is roughly an order of magnitude larger than 
the two other mixing angles and ^23 one can look for models in which 9c is given in terms 
of group theoretical quantities and and ^Is vanish at LO. As shown in Section [6] this can be 
implemented successfully in at least two different ways: a.) one can simply reduce the number of 
Higgs fields in the model by omitting some fields which are allowed to have a non-trivial VEV 
in principle; b.) one can break the dihedral symmetry down to one of its dihedral subgroups, Dq, 
q > 1, instead of Z2. Case a.) has the slight disadvantage that the resulting mass matrices now 
also depend on the choice of the scalar fields and are not only determined by the representations 
under which the fermions transform and the group theory of the dihedral symmetry. Case b.) 
on the other hand does not suffer from this sort of arbitrariness. However, it cannot be realized 
with all dihedral symmetries, since not all of them have dihedral subgroups Dq with q > 1. The 
group D-j which has been used in this paper only has Z2 and Zj as subgroups, since its group 
index is prime. Therefore in the shown examples (for case b.)) the flavor symmetry is taken to 
be Z)i4 instead. The preserved subgroups are of the form D2 =< A^,BA'" >. Also here it is 
necessary to break down to two different D2 groups in the up quark and down quark sector in 
order to generate a non-vanishing Cabibbo angle. One possible choice is m„ = 6 and = 0. 
Finally, we also studied the lepton mixing matrix Vmns numerically. In order to apply the 
results of the mixing matrices found in Section [2] we restricted ourselves to the discussion of 
Dirac neutrinos with a normally ordered spectrum, i.e. the neutrinos have the same properties 
as the other fermions. Since the elements of Vmns are much less constrained than the ones of 
VcKM much more combinations of the group theoretical quantities n, j, mi and rriy can be used 
in order to describe an element of Vmns- However, since we expect that the leptons transform 
under the same flavor symmetry as the quarks, we only considered cases in which the group 
index n is fixed to n = 7 or n = 14. A numerical analysis shows that the experimental fit values 
of the mixing angles can be accommodated well in most of the cases. The strongest constraint 
seems to arise from the upper bound on the reactor mixing angle ^13. Therefore we performed 
fits with two different bounds on siv?{9iz). The results which are shown in Section [7| correspond 
to a very loose bound, sin^(^i3) < 0.1 (which exceeds the 4 a bound [12]), while in the other 
fit the 2a bound, s\v?{9iz) < 0.025, has been used. A common feature of all fits is the fact 
that Jcp vanishes. As shown in Section [7] the condition which minimizes the (13) element of 
Vmns whose absolute value is sin(0i3) also implies Jcp = 0. Furthermore, it turns out that 
the fit parameters 0;, 9i, and a are fixed to a unique value in all cases. In addition to this, we 
also studied how well one could mimic the TBM scenario with a mixing matrix resulting from a 
dihedral flavor symmetry with preserved subgroups. Again, we only considered the cases n = 7 
and n = 14. Our results are similar to the ones found in case of the fit to the experimental best 
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fit values, i.e. a successful fit is possible in several cases. Similar to above, the main restriction 
seems to come from the requirement to pass the bound on sin'^(0i3). For this reason again two 
different bounds on sin^(0i3) have been used. Apart from the cases which lead to similar results 
as above, we observe one additional case, namely if {Vmns)ii is fixed to be cos(^) w 0.7818. 
Unlike in the other cases we can observe CP violation here, i.e. Jcp can take any value between 
—5.776 X 10~^ and 5.776 x 10~^. In contrast to this, the results of the fit of the mixing angles do 
not vary. Unfortunately, sin^(0i3) is very large and therefore a model incorporating this solution 
is disfavored, if contributions from, for example, higher-dimensional operators are not able to 
lower sin^(0i3). 

In the whole discussion we focussed on the case of Dirac neutrinos, since then all formulae found 
in case of the quarks are applicable also to the lepton sector. However, in general neutrinos can be 
Majorana particles. If we assume that they acquire masses from Higgs triplets only, i.e. there are 
no right-handed neutrinos, the possible matrix structures are M5 (Eq.([5])) and a block structure 
(Eq.([3])), see Section [2j Compared to the case of Dirac neutrinos the mixing matrix Ui, is now 
determined by ut U* = diag(mi, m2, "^3) and therefore in general contains Majorana phases. 
This, however, does not matter for the analysis done in Section [71 since there only the absolute 
values of the matrix elements are relevant. Things can change, if we consider the type 1 seesaw 
instead. As we then deal with the Dirac neutrino mass matrix and the right-handed Majorana 
mass matrix, these mass matrices can preserve different subgroups of the flavor symmetry. This 
is, for example, the case in the models [15,16] by Grimus and Lavoura V\ The situation can be 
even more complicated, if the model also includes Higgs triplets. Then all these three matrices, 
i.e. the Majorana mass matrix of the left-handed neutrinos, the one of the right-handed ones and 
the Dirac neutrino mass matrix, can conserve different subgroups of the original symmetry and 
in general no definite statements can be made about the resulting mixing matrix. Furthermore 
we assumed throughout our analysis that the neutrinos have the same mass ordering as all the 
other fermions. However, due to the unknown sign of the atmospheric mass squared difference it 
is also possible that the neutrino mass hierarchy is inverted {m^ < mi < 7712). 
Our study is by no means a complete study of all possible mixing structures which can in principle 
arise from a dihedral flavor symmetry with preserved subgroups. For example, in all cases we 
presented here the subgroups which are preserved in the up and down quark sector have the 
same group structure, i.e. they are either both Z2 or D2 groups. In general, however, these 
group structures could be different. Successful examples which employ subgroups of different 
structures are the A4 (T') models as well as the models by Grimus and Lavoura. As already 
mentioned in the Introduction, in the A4 [T') model the conserved subgroups are Z2 {Z4) and 
in order to predict TBM in the lepton sector. In the first model [15] by Grimus and Lavoura 
the flavor symmetry D4 x ^^""^^ ig broken either to D2, Z2 or is left intact (see [1]). Similarly, 
in their second model [16] with D3 x ^^""^^ as flavor symmetry Z3, Z2 and D3 are the preserved 
subgroups (see also [1]). Both models lead to vanishing ^13 and maximal atmospheric mixing (for 
the leptons). This shows that the usage of subgroups of different group structure leaves much 
more possibilities than the ones shown here. As the complete study of mass matrix structures 
(with det(M) 7^ 0) which can arise from a dihedral symmetry, if a subgroup remains preserved, 
already exists [1], it is only the question how to combine these results in order to get further 
interesting predictions for the mixing patterns of quarks and leptons. One interesting example, 
namely the explanation of the Cabibbo angle, has been studied in detail in this work. 
Finally, let us remark that a common feature of the model(s) shown here and the successful 

^®In their models the Dirac neutrino mass matrix is actually invariant under the whole dihedral symmetry, i.e. 
stems solely from VEVs of fields which transform trivially under the flavor group. 
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(T') models is the need for an additional Z^"""^^^ symmetry which can separate the different 
sectors according to the different subgroups of the flavor symmetry which should be preserved. 
Due to such an additional symmetry an embedding of these models into an 50(10) GUT is in 
general not straightforward. However, assigning the quarks to 

Qx.u\ ~ (li,+l) , ( J ) ' ( I) - (2l,+l) , - , ( I) ~ (^I'-l) (26) 

under x as done in Section [4.1.21 still allows an embedding into SU{5) multiplets. 

Note added: At the final stages of this work the paper [17] by C. S. Lam appeared. He also deals 
with the fact that non-trivial subgroups of some discrete flavor symmetry can help to explain a 
certain mixing pattern and also very briefly mentions that the Cabibbo angle might be the result 
of some dihedral group. 
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A Possible Forms of Vr 



According to the three possible identifications of the eigenvalue c — d there exist three possible 
diagonalization matrices in each sector (up and down sector, charged lepton and neutrino sector) 
U, U' and U" which are shown in Section [3.11 Out of these one can form nine possible mixing 
matrices 1/^^^ = W^" with a, 6 = 1, 2, 3 and Wi G {U, U', U"} where Wi depends on the group 
theoretical phase (pi (the index rrii) and contains the parameters 6i and They all have the 
property that one of their matrix elements, namely the element (ab), is completely determined 
by group theory, i.e. by the index n of the dihedral group, by the index j of the two-dimensional 
representation 2j under which two of three generations of SU (2) l doublets transform and by 
the breaking directions described by mi and m2 in the two different sectors. In the following we 
abbreviate Pi — P2 with a, sm(9i) with Si and cos{9i) with q. 
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The measure of CP- violation Jg!'p is given for the matrices V^\^ as 

Jcp = Jcp{i,(l>i,<l>2;Oi,02,a) , J^p = -Jcp{]AiA2]di,e2,a) , J^p = Jcpij\AiA2\0i,92,a) (27) 
Jc'p = -Jcp{],(l}iA2]di,02,a) , Jc?p ^ JcpH, (1)1, (1)2; 01,02, a) , J'cp ^ -Jcp{],(l>i,(l>2;0i,e2,a) (28) 
J^p = Jcp(j,?!'i,02;6'i,6'2,a) , Jc^p ^ ~Jcp{h(l)i,(l)2;0i,e2,a) , Jg^p = Jcp(j, </'2; 6I1, 6*2, a) (29) 

with Jcp{i,4>i,4'2\di,02,(y) = -g sin((0i - </)2)j) sin(i {(pi - (/)2)j) sin(2 6'i) sin(2 6'2) sin(i {(pi - (p2)i + a) 

(30) 
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Ci 


C2 


C3 
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C5 
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A^ 




A^ 


B 
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2 
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7 




1 


7 


7 
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2 


ll 


1 


1 


1 




1 


1 


l2 


1 


1 


1 




1 


-1 


2l 


2 


2 cos(99) 


2 cos(2(/?) 


2 


cos(3 if) 





22 


2 


2 cos(2(/?) 


2 cos(4(/?) 


2 


cos(6 





23 


2 


2 cos(3 99) 


2 cos(6 99) 


2 


cos(9 ip) 






Table 6: Character table of the group -D7. <^ is 2lL. d are the classes of the group, 
°Ci is the order of the i^^ class, i.e. the number of distinct elements contained in this 
class, °hc; is the order of the elements 5* in the class Ci, i.e. the smallest integer 
(> 0) for which the equation S ^'^^ — 1. holds. Furthermore the table contains one 
representative for each class d given as product of the generators A and B of the 
group. 



B Group Theory of D'^ 

The group D-j has two one- and three two-dimensional irreducible representations which we denote 
as I2, 2-|^, 22 and 23. 1-^ is the trivial representation of the group. All two-dimensional 
representations are faithful. The order of the group is 14. The generator relations for the two 
generators A and B are: 

A^ = 1 , B^ = 1 , ABA = B . 

A and B can be chosen to be 

e~ \ _ / 1 



6 TV i \ 

e~ \ _ / 1 



6 TT i 

e — - 



--ill) 



For the one-dimensional representations 1^ and I2 A and B can be found in the character table 
Table El 

The Kronecker products are: 

X // = /i , I2 X I2 = 1 1 , 12 X 2i = 2j 
[2i X 2i] = ll + 22 , {2i X 2i} = I2 
[22 X 22] = ll + 23 , {22 X 22} = I2 
[23 X 23] = ll + 2i , {23 X 23} = I2 

21 X 22 = 2i + 23 , 2i X 23 = 22 + 23 , 

22 X 23 = 2i + 22 , 
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where /i is any representation of the group and [v x v] denotes the symmetric part of the product 
v X v, while {v x v} is the anti- symmetric one. 

The Clebsch Gordan coefficients are trivial for 1^ x /Lt and I2 x 1.2- 12 ^ ^ non-trivial 
sign appears 



-Ba2 y -1 



for B ^ I2 and ^ ^ ~ 2|. Ij^ and I2 of 2j x 2j are of the form 
ai a2 + a2 a'l '"^ , ai 02 — 02 '-^ I2 



a2 



for ( j , ( '^j 1 ~ 2j. The two-dimensional representations also contained in these products 



read: 



fori = 3: f"2lU2i 
For the rest of the products 2j x 2j we get: 

a2 ) ' ( ^2 ) ^ ( aj62 ) ~ ' ( a2 62 ) ~ ^3 



^1 



All these formulae arc just special cases of the more general formulae given in [1, 10] which hold 
for dihedral groups Dn with an arbitrary index n. 
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C Higgs Potential 



We begin by writing down the potential for the three Higgs fields and -fff 2' which couple 
only to up quarks, i.e. are even under the additional ^^""^^ symmetry. The potential is of the 
same form as V3 in Eq. (ll4p . As mentioned above, it has an accidental U{1) symmetry. 



2 



We have in addition five Higgs fields which are odd under the extra Zg""^'*. These are H^, 
and Xi 2 ■ '^^^ most general potential for these five scalar fields is 

2 



/ 2 \ 2 ,2 \ 2 

. d / rj-dt ^ds 2 . d / \^ rj-dt „d \ , \d I ST^ dt d \ , d / rrd't rrd rjd'f rTd\2 , \di dt d dt d\2 



+ \i\Hf H!^\' +\i\xf xil"" + oi{Hf Hi) (Y^HfHf +at{HfHf) ij^''^^' 



2 \ / 2 

Xi Xi 



i = l 



2 \ / 2 



+ (^E j (^E ^'^j + riiHfHi - HfH})(xfxi - xi'xi) 

+ {ri{Hfxi){Hi\i) + h.c.} + rt (^E \Hfxtt^ + {ri{Hf xi){Hf xt) + h-c} + rt{ \Hfxi\^ + x?l') 

+ {r^{{Hfxt)(xfxi) + {Hfxi)(Xixi)} + h-c} 

+ {ioi{{Hf Ht){Hf xi) + (HfH^mfxf)} + h.c.} + {coUinf Ht){xf H^) + {nf H^){xi^ H^)} + h.c.} 

+ WUiHfxtXnfHi) + {Hfxi){HfHt)} + h.c.} 

This five Higgs potential is free from accidental symmetries. However, the combined potential 
Vu + Vd has an accidental SU (2) x f/(l) x U{1) symmetry. It is broken explicitly by mixing terms, 
which couple the Higgs fields H'^^ 2 ^-i^cl H'^^ 2/Xi 2- 
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The following potential Vmixed contains all such terms, which are invariant under the symmetry 



Knifed = ki{H:^ H:){Hf Hf) + {K2{H:^Htf + h.c.} + Ks\H:^Hff (3 



2 



+ ^,\m^Ht ~ H^^H^f + S[6+9i(i^^r^x?i' + im^xif) + nuAim^ m){Hf Hi) + h.c.} 

{^,,{Hf Hi){Hf Hi) + h.c.} + {kii{Hfxi){Hfxi) + h.c} + ^i2( Ht\^ + \HfH 

fti2( \HVxi? + \Hi\i\') + ^mh:^h:) {Y.HfHf^ + ~ms{H:^H:) (i2xfx 



+ 



if) 



+ Hi 



{ki,{h:^ Hi){H:^ Hi) + h.c.} + {ki^K^xixK^xi) + h.c.} + ^i5( m^nif + m^nif) 
{\H:^xi? + \H:\i\') + ni,{HfHt) {Ylm^Hf^ +{K,,{Hfm){Hfm) + h.c.} + Kis {f^wfrnt 

{t^i^iHfHi) (j^m^^f^ +h.c.} + {n2oiHfHi) (j^Hfm^ +h.c.} 

{k2i{{H:^ Hlt){Hf Hi) + (H:^ Hi){Hf Hi)} + h.c.} + {K22{iH:^ Hi){Hf Hi) + {H:^ Hi){Hf H^)} + h.c.} 

{k23{{h:^ m){Hf Hi) + {h:^ Hi){Hi^ Hi)} + h.c.} + {K24{{H:^Hi){Hi^Hi) + {h:^ Hi){H^^ Hi)} + h.c.} 

{^2^{{Hf Hi){m\i) + {Hfm){Hi^xi)} + h.c.} + {K2e{{Hi^ Hi)ixf Hi) + {nf H^){xf m)} + h.c.} 
+ {^27{{Hfxt){Hi^H^) + {Hfxi){m^Hi)} + h.c.} + {K2s{{H:^Hi){m^xi) + {Hf Hi){Hi\i)} + h.c.} 
+ {K29{iH:^Hi){xfHi) + {H:^ Hi){xf Hi)} + h.c.} + {K;o{{H:^xi){HfHi) + {H:\i)(Hf Hi)} + h.c.} 

+ {^zi{{H:\i)(Hi^ Hi) + (K^ximi^Hi)} + h.c.} + {^32{(//r^^^r)(^^2 V') + iH:^Hi){Hfxi)} + h.c.} 

+ M{HfHi){xfHi) + {HfHi){xf Hi)} + h.c.} 



In our numerical analysis we restricted ourselves to the inclusion of a minimal number of terms 
from Vmixed which break all accidental symmetries such that only three Higgs particles remain 
massless which are eaten by the and boson. As explained in the main part of the text, 
the three terms K2, K5 and K19 are sufficient. 

The numerical example in Section 14.1.11 and Section 14.2.11 needs the following VEV configuration 

{Hf'^) = 61.5 GeV , (Hf) = (H^) = (xf) = (xt) = 61.5 GeV , (H^) = 61.5e-^ GeV 
and {H^) = 61.5 e^ GeV 

which allows real parameters in the potential V^, as all fields Hf, and xf 2 have real VEVs. 
Furthermore we can remove the phase of such that we are left with three complex parameters 
stemming from Vmixed- 

The mass parameters are at the electroweak scale: 

/i^ = 100 GeV , /i^ = 200 GeV , /if = 100 GeV , n% = 200 GeV and /i^ = 150 GeV . 
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One possible setup of parameters is then: 
For Vu we take: 

= 0.959337 , A^* = 2.52548 , A^ = 0.374967 , A^ = -0.588842 , a]' = 1.62353 , 
cr^ = -0.276964 , = -0.283914 . 

For Vd we set: 

Af = 1.70438 , Xf = 3.76598 , Xf = 1.47549 , A^ = -0.344036 , A^ = -0.185157 , A^ = -0.304589 , 

A^ = -0.733236 , af = 0.22429 , af = 4.6792 , = -0.87457 , = -2.0284 , = 0.961454 , 

ai = 0.649984 , rf = 2.96557 , = 1.22903 , ri = -2.02133 , = -1.22242 , = -2.31577 , 
r| = 2.38236 , = -0.660102 , uf = 0.452165 , = -2.112 , = -1.63452 . 

and for the three complex nonzero couplings from Vmixed- 

K2 = -0.638073 + i 0.0277608 , K5 = 0.312782 + i 0.140162 , K19 = -0.278402 - i 0.124756 . 

Note that all parameters have absolute values smaller than 5 and hence they are still in the 

perturbative regime. 

With these parameter values we obtain the desired VEV structure. 
The Higgs masses are then 

513 GeV, 499 GeV, 426 GeV, 414 GeV, 386 GeV, 365 GeV, 321 GeV, 266 GeV, 246 GeV, 227 GeV, 
178 GeV, 159 GeV, 134 GeV, 81 GeV and 55 GeV 

for the neutral scalars. Due to the explicit CP violation in the potential we can no longer 
distinguish between scalars and pseudo-scalars. For the charged scalar fields we get 

367 GeV, 333 GeV, 294 GeV, 261 GeV, 145 GeV, 115 GeV and 55 GeV. 

They are therefore in general too light to pass the constraints coming from direct searches as 
well as from bounds on FCNCs. Nevertheless, soft breaking terms of mass dimension two of the 
order of 10 TeV could lift the masses above these experimental bounds. 
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